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Abstract 

A family of boundary conditions corresponding to exclusion processes is 
introduced. This family is a generalization of the boundary conditions 
corresponding to the simple exclusion process, the drop-push model, and 
the one-parameter solvable family of pushing processes with certain rates 
on the continuum [1-3] . The conditional probabilities are calculated using 
the Bethe ansatz, and it is shown that at large times they behave like the 
corresponding conditional probabilities of the family of diffusion-pushing 
processes introduced in [1-3]. 



PACS numbers: 05.40.-a, 02.50.Ga 



1 Introduction 



In recent years, the asymmetric exclusion process and the problems related to 
it, including for example the bipolimcrization [4], dynamical models of interface 
growth [5], traffic models [6], the noisy Burgers equation [7], and the study of 
shocks [8,9], have been extensively studied. The dynamical properties of this 
model have been studied in [9, 11]. As the results obtained by approaches like 
mean field are not reliable in one dimension, it is useful to introduce solvable 
models and analytic methods to extract exact physical results. Among these 
mcthodes is the coordinate Bethe-ansatz, which was used in [12] to solve the 
asymmetric simple exclusion process on a one-dimensional lattice. In [1], a sim- 
ilar technique was used to solve the drop-push model [13], and a generalized 
one-parameter model interpolating between the asymmetric simple exclusion 
model and the drop- push model. In [2], this family was further generalized to 
a family of processes with arbitrary left- and right- diffusion rates. All of these 
models were lattice models. Finally, the behaviour of latter model on continuum 
was investigated in [3] . The continuum models of this kind are also investigated 
in [14, 15]. In [16], another generalization of boundary conditions correspond- 
ing to these processes was introduced, which shows an exclusion process with 
annihilation as well as diffusion. In [17], the Bethe-ansatz method was used 
to investigate reaction-diffusion processes involving particles of more than one 
type, on a lattice. 

In the generalized model interpolating between the asymmetric simple ex- 
clusion model and the drop-push model [1-3], there are two parameters A and fi, 
which control the pushing rate. Normalizing the diffusion rate to one, it is seen 
that the sum of these two parameters should be one to ensure the conservation 
of probability. These two parameters appear in the boundary condition used 
instead of the reaction. The question is that on continuum, what other kind of 
boundary conditions, corresponding to other reactions, can be imposed. This is 
what is investigated in the present paper. 

The scheme of the paper is the following. In section 2, the allowed boundary 
conditions are investigated. The criterion to choose a boundary condition is 
the conservation of probability. It is shown that the boundary condition is 
characterized by a one-variable function which vanishes at the origin. 

In section 3, the Bethe-ansatz solution for the 2-particle probability of this 
process is obtained, and its large-time behavior is investigated. It is shown that 
the large-time behavior of the system is determined by only the first nonzero 
derivative of the function determining the boundary condition, at the origin. If 
this function does have a nonvanishing linear term, then that term determines 
the large-time behavior of the system, and the results of [3] are recovered. 

Finally, in section 4, boundary conditions are investigated for which the first 
derivative of the function determining the boundary condition vanishes at the 
origin. The mean position of the particles, and the effective diffusion parameter 
are obtained. 
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2 Boundary conditions 



Consider a collection of N particles diffusing on a one-dimensional continuum. 
So long as the particles do not encounter each other, the master equation gov- 
erning the probability density of finding these at xi < • • • < xn is [3] 

|p(x;t) = lv 2 P(x;i). (1) 

In [3], it was shown that the exclusion pushing condition can be written as a 
boundary condition with two parameters (in fact with only one independent 
parameter, since A + /i = 1): 



^ dxi ^ dx i+ i ^ 



= 0. (2) 



To obtain this system, in [3] it was started with a collection of particles diffusing 
on a one-dimensional lattice. Each particle diffuses to the right with the rate 
one, if the right site is free. If the right site is occupied, the particle may push 
other particles and go to the right, with a special rate. That is, we have the 
following process. 



AA---A$^$AA---A, with the rate r n , (3) 

where 

r " := [ i+ '" + (?)T- (4) 

It was shown in [1], that such a system can be described by the master equation 

P(x;f) =^P(x-e i ;t)-iVP(x;i) ) a:i < • • • < x N , (5) 

i 

subject to the boundary condition 



P(...,X{ = x, Xi+i = x, . . . ; t) =AP(. . . ,x % = x, x l+1 = x + l,...;t) 

+(iP(. ..,Xi=x-l, x i+x =x,...;t). (6) 

In (JjjJ, ej is a vector the components of which are all zero, except for the z'th 
component, which is one. In [3], and J5J were obtained as the limit of (0 
and ((HI when the probability is a slowly- varying function, with a suitable Galileo 
transformation, so that one has diffusion to right and left (with equal rates). So, 
JQl and J5J) describe a system of particles diffusing on a continuum and pushing 
each other with certain rates, as introduced in J2J and (0J. 
For the special case N = 2, (JIJ and J2J become 

^P{x u x 2 ;t) = ^{dl + dl)P{x u x 2 ;t), (7) 
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and 

fxdiP(x, x) = \d 2 P{x, x). (8) 

Using 

dP{x,x) =d 1 P(x,x) + d 2 P(x,x), (9) 
it is seen that the boundary condition JSJ is equivalent to 

{di-d 2 )P(x,x) = {\-(j.)dP(x,x). (10) 

The passage from (j5J) to (J2J) is not unique. But any boundary condition has 
to ensure probability conservation. Consider the evolution (JJJ, and note that it 
is valid only for the physical region X\ < x%. Integrating this over the physical 
region, one arrives at 



d 
~dt 



J dxidx 2 P(xx, x 2 ) = ^ J dx(di - d 2 )P(x,x). (11) 

But using i|l(J|) . it is seen that the right-hand side of (|11|) is zero. So, 

d f 

— I d Xl dx 2 P(x u x 2 ) = 0. (12) 

This is the conservation of probability. However, the only thing needed to ensure 
this conservation is that (pi — d 2 )P(x, x) be a total derivative of some function 
with respect to x. So, one can in general use a boundary condition like 

(d\ — d 2 )P(x, x) — a\dP(x, x) + a 2 ddP(x, x) + ■ ■ ■ 

= f(d)P(x, x), (13) 

where / is an analytic function satisfying 

/(0) = 0. (14) 

Such a boundary condition ensures probability conservation. In [16], a slightly 
different boundary condition was introduced, which in this language is equiva- 
lent to /(0) < 0. It was seen there such a boundary condition describes annihi- 
lation as well as diffusion. 

3 The Bethe-ansatz solution 

The Bethe-ansatz solution to the Master equation QJ, with the boundary con- 
dition l|13|) is 

P(x;t) =e £t *(x), (15) 

where "J satisfies 

B*(x) - ^V 2 v]/(x), (16) 



3 



and the boundary condition l|13fl . The solution to these, is 

* k (x) =J2A*e^>x, (17) 

a 

with 



\Y. k l ( 18 ) 



E ~ 2 

i 

provided one can find A^s so that (|13fl is satisfied. The summation in l|17fl is 
over the permutations of N objects, Applying (| 1 31) to the case Xj — Xj+i, one 
obtains 

i[k a(j) - k a(j+1) ]{A a - A aaj ) = i[f(K(j)) + f(K(j+i))]( A a + A a<Xj ), (19) 
where / is defined as 

f(k) := -if(ik), (20) 

and Uj is that permutation which changes j to j + 1 and vice versa, and leaves 
other numbers unchanged. From (|19|l . A a can be obtained as 

Aaa-j = S(k rT (j),k a {j +1 ))A (7 , (21) 

where 

s(k 1 ,k 2 ) = k2 - h + Kkl + k2) 

ki - k\ - f(ki + k 2 ) 

= 1 + 2/(fcl + fc2) . (22) 
ki - ki - f(h + k 2 ) 

As cria2<Ji = cr 2 aicr 2 , the coefficients A, should satisfy the consistency condition 

A(j 1 a 2 &i — A ( j 2( j 1 (j 2 ■ (^^) 

The above criterion, in terms of S is 

S(k 2 , k 3 )S(k 1 ,k 3 )S(k 1 ,k 2 ) = S(k 1 ,k 2 )S{k 1 ,k 3 )S(k 2 , k 3 ), (24) 

which is obviously an identity. 

Using this Bethe-ansatz solution, the conditional probability can be written 

as 

P(x; t|y; 0) = J ^-p tt k (x) e E ^-^, (25) 

where in 'I'k, the coefficient of e lk x is set to be equal to one. 

This shows the integrability of the system, in the sense that the iV-particle 
scattering matrix (A^'s) can be expressed as a product of two-particle scattering 
matrices. 



4 



3.1 The conditional probability for the 2-particle sector 

The conditional probability for the 2-particle sector, is written as 

P(x;t|y;0) = J dhdk 2 [e lklXl+lk2X2 + S(k u k 2 )e lklX2+lk2Xl ] e Et ~ lk y 

(26) 

so that it satisfies (0 and the the initial condition 

P(x,0\y,0)=6(x 1 -y 1 )6(x 2 -y 2 ) (27) 

for the physical region (xi < x 2 , yi < y 2 ). 

Using the form of S, one can write this conditional probability as 



P(x; t\y; 0) =J 0e Bt -* k ' y \j( k ^+ k ^) 



2f(h + k 2 ) 



k 2 - ki - f(h + k 2 ) 



e i(kiX2+k2Xi) I (28) 



In the above integration, the possible ambiguity arising from the the pole of the 
fraction is removed through k\ — ► k\ + i e and k 2 — > k 2 — i e, where the e — ► + 
limit of the integral is meant. This ensures that the probability tends to zero 
as x\ — ► — oo or x 2 — ► +oo. 

The two first integrals of the right-hand side are easily calculated. So, 

P(x;i|y;0) = — { e -[(^-w) 2 +(^-^) 2 ]/(2t) + e ~(z 2 +z 2 )/(2t) j + ^ (2Q) 
where 

z 1 :=x 1 -y 2 , z 2 :=x 2 -yi. (30) 
To obtain the third integral (I3), one uses the change of variable 

k := fci + k 2 , q := k 2 — hi — f{k\ + k 2 ). (31) 



Then, 



where 



7 3 = / dMq c -ttk 2 + \ Q +Mk)] 2 }/4 c itkA 2 + \ Q +J(k)]A 1 }/2 

47T 2 q 

M dAe -[A 2 +2tAfm]/(4t)^ (32) 

-00 

A\\=z\ — z 2 , A 2 :=zi + z 2 . (33) 



5 



Simplifying \.'>'2\) . one arrives at 



tk 2 - 2ikA 2 A 2 
x e 4 e ~ 4< 

Ax — A ( d \ A 2 + A 2 

To show that 1291) does indeed satisfy the initial conditions, first note that 

P(x; 0|y; 0) = 5{ Xl - y x ) S(x 2 - y 2 ) + S(xx - y 2 ) 5{x 2 - yi) + h{t = 0). (35) 

The second term is obviously zero in the physical region x\ < x 2 and y± < y 2 . 
For the third term, we have from the first equality in i|32|) 

U t = 0) = I % ](k) e ^+M](Km f Jl_ (36) 
J 47H J q-ie 

The integral over q is, however, zero in the physical region. Since in the physical 
region A\ < 0, and one can close the integration contour by adding a large 
semicircle in the lower half plane of q, to the real line. There are no singularities 
inside the closed contour, so the integral vanishes. So Q29JI does indeed satisfy 
the initial conditions. 

Now consider the large-time behavior of the system, using the change of 
variable 

Ai =: 2a,Vt, (37) 

the integral I3 becomes 

h _ 1 r da vis (1.1.) f t{ai - a)l e -<.- + .i». (38) 

It is seen that at large times, the dominant term comes from the first nonzero 
term in the expansion of /. This means that if the first derivative of / does not 
vanish at the origin, then at large times, / is equivalent to a linear function. 
But a linear / is just the boundary condition used in [3]. So at large times, 
the boundary condition has effectively only one free parameter determining the 
interaction. 

If / is an at-most-quadratic polynomial, then the integration over k in 
can be done. For the simple example 

f(x) = a 2 x 2 , (39) 
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the result would be 

<- Ai dA r a 2 
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_ DO 2ir \ [t + 2a(A - Ax)} 2 t + 2a(A - Ax) 

| ^2 

e _ tt _ 31+8^17^77 _ (40) 



y/t[t + 2a{A- Ax)} 

4 Boundary conditions corresponding to func- 
tions / with vanishing first derivative at the 
origin 

Taking the form 

f(x) = a n x n + 0(x n+1 ) (41) 

for / (where n > 1), it is seen from l|38|l that at large times, the leading term 
of I 3 is 

f a i fin 

d*-Le- (a2+4) - (42) 



7iVt n+1 J-oo da 2 l 

The integration in the left-hand side, can now be easily done, and one arrives 
at 

7 ^?visr [1 + erf(< " )l S^- (43) 

Two important quantities to be calculated, are the mean position of the 
particles, and the variance of the mean position. The expectation value of any 
function g(xx,x 2 ), is obtained through 

(g(xx,x 2 )) := d 2 x P(x.\y)g(xx,x 2 ), 

J Xl<X2 

=(9( x i> x 2))o + (g(xx,x 2 )) 3 , (44) 

where 

(g(xx,x 2 ))3 ■= / d 2 x I 3 (xx,x 2 )g(xx,x 2 ), 

/oo /•(»l-Va)/(2Vi) 
da 2 dax h{xx,x 2 )g(xx,x 2 ), (45) 

-oo J — OO 

and 

{g(xx,x 2 )) Q := f d 2 x -L e -^-y^+^-^^g(xx,x 2 ), 

Jx 1 <x 2 ZlTt 

+ [ d 2 x — e-^-y^ 2+ ^--y^lWg{x 2 ,Xx). (46) 
Jx 1 >x 2 27rt 



7 



If g(x\ 1 X2) — g(x2,xi), then the above expression becomes simpler: 

(g(x 1 ,x 2 )) = / d 2 x le-K^^^l/Mj^, X2 ). (47) 
Using these, we want to calculate 

X:=±(xi+x 2 ), (48) 

and 

(AXf:=~(x 2 +x 2 2 )-X 2 . (49) 
Using l|47[l. it is easily seen that 

^ yi +2/2 ,_ n s 
*o = g , (50) 

and 

i(x?+a ; |)o = (Xo) 2 +t. (51) 

From (|38|l . it is seen that ^3 is the n'th derivative of some function with 
respect to a 2 . So, if g is a polynomial of order less than n with respect to a 2 , 
the right-hand side of l|45|l vanishes; this can be seen through n times integration 
by parts. The functions relevant to the calculation of the mean position and 
the variance, in terms of a\ and a 2 , are polynomials of a\ and a 2 of at most the 
second degree. From (I41|) . it is seen then that only (03)3 may be nonvanishing, 
and that even this term vanishes if n > 2. So, for n > 2 we have X = X and 
Var = Varo. For n — 2, it is seen that I3 contains second and higher derivatives 
of functions with respect to a 2 . In calculating (a|) 3, only the second derivatives 
are relevant. That is, one can write 

h = -^[l + erf( 0l )]|^e-°2 +h, (52) 
where I3 is irrelevant to the calculation of (a?,) 3. So, one arrives at 



2a 2 

\ a 2)3 = 

Using 



(»l-l/a)/(2Vt) 

oi> 3 = / rfax [1 + erf (5+) 



1? + s| = + (54) 



where -B is a first order polynomial, one arrives at 
~(xl + xl) 3 =t(al) 3 



f(Vl-V2)/(2Vt) 

2a 2 Vi / da x . (55) 
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From this, l|5ip. and the fact the mean does not depend on / (so long as the 
first derivative of / vanishes at the origin) it is seen that 

r(vi-V2)/(2Vt) 

Var = t + 2a 2 Vt da x [1 + erf(ai)]. (56) 

This becomes even simpler at large times. At large times, the upper limit in 
the above integral tends to zero. So, using 

1 

[l + ei{(x)]dx = -=, (57) 
. -OO V 7 *" 

one arrives at 

V&T = t + 2a 2 ^/tj^+0{r 1/2 ). (58) 
The effective diffusion parameter is defined as the derivative of the above: 

- d r = 1 + vk + oit f ) - (59) 

It is seen that as t — * oo, this parameter tends to one, in agreement with [3, 12]. 
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